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ABSTRACT 

The origin, evolution and role of magnetic fields in the production and 
shaping of proto-planetary and planetary nebulae (PPNe, PNe) is a subject of 
active research. Most PNe and PPNe are axisymmetric with many exhibiting 
highly coUimated outflows, however, it is important to understand whether 
such structures can be generated by isolated stars or require the presence of a 
binary companion. Toward this end we study a dynamical, large-scale a — Q 
interface dynamo operating in a 3.0 Mq Asymptotic Giant Branch star (AGE) 
in both an isolated setting and one in which a low-mass companion is embedded 
inside the envelope. The back reaction of the fields on the shear is included 
and differential rotation and rotation deplete via turbulent dissipation and 
Poynting flux. For the isolated star, the shear must be resupplied in order 
to sufficiently sustain the dynamo. Furthermore, we investigate the energy 
requirements that convection must satisfy to accomplish this by analogy to 
the sun. For the common envelope robust dynamo results, unbinding 

the envelope under a range of conditions. Two qualitatively different types of 
explosion may arise: (i) magnetically induced, possibly resulting in collimated 
bipolar outflows and (ii) thermally induced from turbulent dissipation, possibly 
resulting in quasi-spherical outflows. A range of models is presented for a variety 
of companion masses. 

Subject headings: planetary nebulae: general - stars: AGB and post-AGB - 
stars: low-mass, brown dwarfs - stars: magnetic fields 

1. Introduction 

Most planetary and proto-planetary nebulae (PNe, PPNe) are highly aspherical 
and exhibit a diversified morphology of axisymmetric structures and/or collimated jets. 
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However, the production of such richly varied systems and shaping mechanisms, remains 
an open topic (Bahck & Frank (2002)). A central question is whether a binary is required 
to produce such asymmetries or if an isolated Asymptotic Giant Branch star (AGE) is 
sufficient. 

Recently, detection of magnetic fields in AGB stars (Etoka & Diamond (2004), Bains 
et al. (2004)) and the central stars of PNe (Jordan et al. (2005)) has sustained interest in 
magnetic launching and collimation mechanisms. Observational evidence of a magnetically 
collimatcd jet in an evolved AGB star (Vlemmings et al. (2006)) further suggests that the 
magnetic field may play a dynamical role. 

Single star magnetic outflow models have been proposed as mechanisms for producing 
and shaping PPNe and PNe (Pascoli (1997), Blackman et al. (2001)). Whether such 
models can power and shape the PPNe is uncertain (Soker (2006)). In particular, envelope 
dynamos are expected to be short (< 100 yrs) and drain differential rotation energy 
rapidly (Blackman (2004)), making it unlikely that isolated stars can produce observed 
asymmetries. If convection can resupply differential rotation energy, then an envelope 
dynamo in an isolated AGB star may be viable. Anisotropic convection in the sun resupplies 
differential rotation through the A-effect (Riidiger (1989), Riidiger & HoUerbach (2004)), 
however it is not clear if such a mechanism operates in AGB stars. 

Rather than single star models, the observed asymmetries may instead be the result 
of energy and angular momentum supplied by binary interactions. This is supported by 
recent surveys suggesting that most, if not all PNe involve binary systems (De Marco et 
al. (2004), Sorensen & PoUacco (2004), Moe & De Marco (2006), Mauron & Huggins 
(2006)). Although there are many different types of binary interactions and outcomes, 
here we focus on common envelope (CE) evolution in the context of PNe progenitors 
(Ibcn & Livio (1993)). A common envelope (CE) model in which low-mass (< 0.3 Mq) 
companions were embedded into the envelope of a 3.0 Mq AGB star was investigated in 
Nordhaus & Blackman (2006). The CE evolution is advantageous as it can supply angular 
momentum in an extremely short period (< 1 yr) and produce a range of PPNe outflows. 
Such a model would predict white dwarf + brown dwarf close binaries that survive the 
CE phase. Recently, a WD+ BD binary in a close orbit has been detected (Maxted et al. 
(2006), Burleigh et al. (2006)). A separation distance of 0.65 Rq indicates that the system 
incurred a common envelope phase in which the brown dwarf was responsible for ejecting 
the progenitor envelope. This further motivates more detailed study of CE induced PNe. 

In this paper, we reinvestigate the magnetic model presented in Blackman et al. (2001) 
in more detail, in an effort to determine the viability of a single star dynamo. We compare 
the results to a model in which the rotation profile is supplied by a CE phase as in Nordhaus 
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& Blackman (2006). In section 2, we review previous work, compare single star evolution 
to that resulting from a CE interaction and focus on the depth that the poloidal field can 
diffuse into the shear zone. We present our interface dynamo in section 3, including the 
detailed back reaction of the fields on the shear, the generation of heat through turbulent 
dissipation and the spin down of the star due to Poynting flux. In Section 4, we present the 
results of our model for three cases: (i) an isolated dynamo in which convection does not 
resupply differential rotation, (ii) an isolated dynamo in which a fraction of the convective 
energy resupplies differential rotation energy and (iii) a dynamo resulting from the in-spiral 
of a low-mass (< 0.05 Mq) companion inside the stellar envelope. We conclude in section 5. 



2. Dynamos, Common Envelopes and Isolated AGB Evolution 

A central issue in magnetic PNe progenitor models is whether an isolated AGB star 
can sustain the necessary field strengths and corresponding Poynting flux to unbind the 
stellar envelope and produce coUimated outflows. Several authors have appealed to various 
mechanisms with which to produce magnetically mediated outflows in isolated settings 
(Tout & Pringle (1992), Pascoh (1993), Pascoh (1997), Blackman et al. (2001), Soker & 
Zoabi (2002), Matt et al. (2004)). Soker & Zoabi (2002) appeal to em a'^ - Q dynamo 
operating in the AGB envelope as a means of enhancing dust formation near magnetic cool 
spots on the stellar surface. The corresponding fields arc not strong enough to dynamically 
alter the geometry, however enhanced mass-loss near the spots could form an elliptical PN. 
Such a model docs not produce strongly bipolar PNe. 

Other authors have investigated a — Q dynamo models which tap into the differential 
rotation energy reservoir between core and convective zone. Pascoh (1997) solved for a 
steady-state, radial dynamo model from inside the AGB core to produce fields throughout 
the envelope. Toroidal field strengths of 10^ G are obtained at the surface of the core with 
poloidal field strengths about an order of magnitude lower. No back-reaction of the fields 
on the rotation profile was included. 

Blackman et al. (2001) investigated a simplistic interface dynamo model (Parker 
(1993), Thomas et al. (1995)) operating at the base of the convective zone in our 3 Mq AGB 
star (see Fig. 1). Angular momentum is conserved on spherical shells as the star evolves 
off the main sequence and the resulting rotation profile is used to calculate field strengths. 
To drive PPNe, the corresponding dynamo must operate through the entire lifetime of the 
AGB phase (~ 10^ yrs) until radiation pressure has bled most of the envelope material 
away. Only then can the Poynting flux unbind the remaining material. But, there are 
challenges for this model. The differential rotation zone is chosen to be ~ 1/2 of the total 
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distance from core to convective zone and ~ 1/2 the length of the convective region. This 
results in only a small fraction of the free shear energy available for field amplification. The 
majority of the shear energy is located deeper, near the core-envelope boundary. The actual 
thickness of the differential rotation layer tapped by the dynamo is uncertain. If more of the 
shear energy were extracted, the envelope could be blown off too early. On the other hand, 
any dynamo operating beneath the envelope would drain differential rotation on time scales 
short compared to the AGB hfetime. Only if the differential rotation is re-seeded might 
this problem be overcome. Even in the sun, the complex interaction between anisotropic 
convection and the resupply of differential rotation is not fully understood. Recently, 
the role of downward pumping and penetration depth in the solar tachocline has been 
investigated (Browning et al. (2006), Dikpati et al. (2006)). 

If isolated star models fail to generate sufficient Poynting flux, common envelope 
evolution provides an alternative mechanism with which to supply significant differential 
rotation energy over very short periods (~ 1 — 10 yrs). The in-spiral of a low-mass secondary 
(< 0.3 Mq) through the envelope of a 3.0 Mq star in the AGB phase was investigated in 
Nordhaus & Blackman (2006). Three qualitative scenarios were found dependent on the 
mass of the companion: (i) direction ejection of envelope material resulting in an equatorial 
outflow, (ii) spin-up of the envelope resulting in an explosive dynamo driven jet along the 
rotational axis and (iii) tidal shredding into a disc which facilitates a jet. In this paper, we 
investigate (ii) further, in addition to presenting results for an isolated star dynamo. 



In order to determine the temporal behavior of the large-scale magnetic field, we 
employ the mean-field induction equation which results from averaging the standard 
induction equation in the presence of hehcal velocity fluctuations. The result is (Parker 



where U is the mean velocity field, B the mean magnetic field, A the micro-physical 
magnetic diffusivity, f3 the turbulent diffusion such that X <^ f3, a the pseudo-scalar helicity 
coefficient, and S =< u x b aB -|- /?V x B the turbulent electromotive force (Moffatt 



3. Dynamical Equations 



(1979)) 




(1) 



(1978)). 
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Although we envision the dynamo engine operating in a spherical or quasi-spherical 
AGB star, for present purposes we work in local Cartesian coordinates. Such interface 
dynamo models have been employed in a variety of systems ranging from late type 
stars (Robinson & Durney (1982)) to white dwarfs (Markiel et al. (1994)) to supernova 
progenitors (Blackman, Nordhaus & Thomas (2006)). The coordinate system and global 
geometry is presented in Fig. 1. The convection zone extends from the stellar surface to the 
interface at r — rc- In this layer, convective twisting motions convert buoyant toroidal field 
into poloidal field through the ct-effect. Below the convection layer, the differential rotation 
zone shears poloidal field back into toroidal field through the Q-effect. Defining the vector 
potential as A = {Ax,A,Az) and decomposing the mean field into toroidal and poloidal 
components, B = {0, B,dxA), generates two coupled equations for the time evolution of 
both components of the magnetic field. 

In order to capture aspects of the 2-D geometry within the framework of our simple 
Cartesian 0.5-D model, we break the turbulent diffusion into two distinct values: 
corresponding to the poloidal field which grows primarily in the convective region and 
(3^, corresponding to the toroidal field which is amplified in the differential rotation zone. 
We also employ (3^ as the turbulent diffusion coefficient for the toroidal velocity. Since 
the convective region is highly turbulent and the differential rotation zone is more weakly 
turbulent, we have that /^^ <C Pp. The ratio of these can be defined as the turbulent 
magnetic Prandtl number as Pvp = Then, assuming axisymmetry (i.e. dyS — 0) for all 
mean quantities S and defining the velocity field as U = (0, U, u), we obtain 
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d^ad^A + d^Ad.U - ud,B - uB/L + (pi + dfj B (2) 
dtA^aB + Pp(dl + df)A. (3) 

where Bd^u ~ uB/L represents a buoyant loss of magnetic fiux and u > Q {to be made 
exphcit later) . We further assume that the Fourier transforms of the fields are proportional 

to delta functions implying that the mean-field has one large scale. We correspondingly 



define 
time. ' 



B,A 



[B{t),A{t)]e 



i(kxX+k, 



'■^^ where A{t) and B{t) are complex valued functions of 



hen, setting k'^ — kl + and using 5^C/ ~ —rcAfl/L yields the following 



d,B ^ak^A- ik^r^A— - iuk.B - uB/L - p,k B 

dtA = aB- /Spk'^A, 



(4) 

(5) 



where the rotation profile across the differential rotation layer varies from Q at the interface 
to + AQ at Tc — L. Thus, AQ is a measure of the shear in the differential rotation zone. 
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If AQ = 0, then the system exhibits sohd body rotation. In addition, we parameterize 
the turbulent diffusion coefficients as [/3</,,/3p] = [c<^,Cp]vLi. where and Cp are distinct 
dimensionless constants and v is a typical convective velocity in the a-layer. 

For the loss of toroidal flux due to magnetic buoyancy, we use the following expression 
for the rise velocity of a magnetic flux tube (Parker (1955), Thomas et al. (1995)): 

3g /a\2 2>QV^^ 



^ 8 VlJ 47rpv 32 v 

where a (assumed to be L/2) is the radius of the flux tube, V a is the Alfven velocity 
associated with the large scale field and Q is a dimensionless constant of order unity. 



3.1. Evolution of Vt and AO 

As differential rotation energy is tapped by field amplification, the corresponding 
Poynting fiux drains rotational energy at the interface. In addition, turbulent diffusion 
converts differential rotation energy into heat which may also be used to unbind the AGB 
envelope. Therefore, to investigate the interaction between field amplification, differential 
rotation and rotation, we derive equations for the evolution of AQ and Q,. The mean-field 
Navier-Stokes equation is given by 

a^u^-u- vu + ^(b- v)b + /3^v^U (7) 

where p is the fluid density, /J^ the turbulent viscosity and we have assumed that b <^ B 
and u ^ U in the weakly turbulent shear layer. Then, taking the ^-component (see Fig. 1) 
yields the following 

dtU ~ -^d^Ad.B + 13^ (dl + dl) U. (8) 

Using the fact that, dzU — rcAQ/L, we can hnk AQ with U as follows 

-dt {rA^) = dt (y (re) -U{rc- L)) ~ dt [idjl) . (9) 

Subtracting the time-dependent velocity equation at Tc — L from Tc and using the relation 
dtU\r^ — dtU\r^-L ~ LdzdtU\r^ yields the following 
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dtAn = — ^ k^kl (Re(A)Re{iB) + Re{iA)Re{B)) - ^k^kMiA)Re{iB) - ^Ail, 

(10) 

where we have assumed p — p and thus d^p ~ {p2 — Pi) / L is the change in density across 
the shear layer. 

In addition to the dynamic shear, we allow for the rotational energy of the field- anchored 
matter to drain via Poynting flux. No appreciable toroidal field amplification occurs above 
the interface, so we calculate the Poynting flux at the base of the convection zone. The 
total integrated Poynting flux at Tc is given by 



Lmag = ^ 1 (E X B) ■ dSe 

~ -Re(B)Re{d^A)nrl. (11) 

When the toroidal and poloidal fields are out of phase, there is a maximum magnetic 
luminosity that is not the respective product of the maximum individual field strengths 
(Blackman, Nordhaus & Thomas (2006)). This is a generic feature of our interface dynamo 
solutions and will be presented in Section 4. 

To arrive at a dynamical equation for Q, we must calculate the available rotational 
energy in the shear layer. We estimate the rotational energy in the differential rotation 
zone as Erot ~ M/^qt'^Q'^ /2. Various mass secondaries can supply a range of differential 
rotation in addition to rotational energy. For a 0.05 Mq brown dwarf, the total shear 
energy from the convective boundary to the tidal shredding radius is ~ 5 x 10"^^ ergs in 
the AGB envelope. This is approximately 4 times the binding energy of the entire AGB 
envelope. As the star evolves into its thermal pulsing phase, the convective zone deepens 
and both the thickness and mass of the shear layer shrink. It may also be possible to power 
an interface dynamo during a later phase in the stars evolution. In this paper, we focus on 
the beginning of AGB phase. 

Equating the time derivative of the rotational energy with the magnetic luminosity 
gives 



a.n MB)Re(a.A)Lr. 

In arriving at Eq. (12), we have multiplied the available rotational energy by a factor of 
6/L. The penetration length, 5, represents the depth at which the poloidal field can diffuse 
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into the shear layer. If 5/L — 1, the total shear energy in the differential rotation layer 
is available for extraction. We estimate the penetration depth as the distance that the 

poloidal field can diffuse into the shear layer during a cycle period. In the kinematic limit, 

1 

the cycle period is given as r = 27r ( agAriofcr ) ^ '^here and AQq are the initial values 
(Blackman, Nordhaus & Thomas (2006)). The cycle period does increase in the dynamical 
regime, however it does not change appreciably from its initial value. Therefore, r serves as 
a lower limit for the cycle period. Correspondingly, we define the penetration depth as 

(13) 



3.2. Evolution of a 



In addition to a dynamical equation for shear, a quenching can be understood through 
magnetic helicity conservation (Kleeorin & Ruzmaikin (1982), Blackman & Field (2002), 
Brandenburg & Subramanian (2005)). In the absence of boundary terms, magnetic helicity 
is well conserved and the build up of A ■ B corresponds to a build up of large-scale field. The 
small-scale helicity then grows to equal and opposite magnitude of the large-scale helicity. 
To maintain simplicity, here we appeal to a parametrized form of a which approximates 
the non-hnear quenching (Blackman, Nordhaus & Thomas (2006)). We adopt the following 
profile 



a — aQ (n/Qo) Exp 
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B I^Ti 



(14) 



where ao = Cq,-^^^, Cq, a dimensionless constant, Vlo the initial rotation rate at the interface 
and pi the density in the middle of the convective region. 



4. Numerical Results 



To investigate various interface dynamo configurations, Eqs. (4), (5), (10) and (12) are 
solved numerically, the solutions of which represent the time evolution of B^^ Bp, AQ and 
Q. In each case, we employ a 1 G seed field for the real components of both the toroidal 
and poloidal field. The fields grow until they are quenched through a drain of the available 
differential rotation energy. In all cases, the saturated dynamo is r2-quenching dominated 
and not a-quenching limited in contrast to Blackman et al. (2001). 
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We focus on three types of shell dynamos: (i) that of an isolated AGB star, (ii) that 
of an isolated AGB star in which convection resupplics differential rotation, (iii) that of 
an AGB star which has been spun up by a companion in-spiraling through the stellar 
envelope. We use independent radial rotation profiles for the above cases to calculate the 
initial value of AQ. For the isolated AGB star, we assume that as the star evolves off 
the main sequence, angular momentum is conserved on spherical shells yielding a rotation 
profile (X r~^. We consider this case in detail in Section 4.1. Because convective energy 
may resupply differential rotation energy analogous to what occurs in the sun, we consider 
a range of resupply rates in Section 4.2. In Section 4.3, we consider a rotation profile 
generated from the in-spiral of a low-mass secondary through the stellar envelope. The 
in-spiral time is short compared to the AGB lifetime and rapidly creates a strong shear 
region beneath the convective zone. In addition, the in-spiral time is less than or equal to a 
cycle period, r, so the angular momentum and energy are supplied to the dynamo on time 
scales short compared to its growth time. 

4.1. Isolated Dynamo Without Reseeding AQ 

It has been suggested (Blackman ct al. (2001)), that isolated AGB stars may be 
able, through dynamo action, to unbind their envelopes and produce collimatcd outflows. 
As the star expands onto and through the AGB branch, the core contracts while the 
envelope expands, creating a shearing profile throughout the interior. A differential rotation 
zone, coupled with the above convective region provides conditions for large-scale field 
amplification. However, field growth requires draining differential rotation energy. In the 
sun, the rotation profile is re-established through a transfer of convective energy by way of 
the A-effect (Riidiger & HoUerbach (2004)). This results in a steady-state rotation profile 
and a magnetic cycle with quasi-steady peak field. However, if convection does not resupply 
differential rotation energy, any resultant dynamo would be a transient phenomena. 

Blackman et al. (2001) investigated an isolated shell dynamo operating in the AGB 
phase of our model star. The depth of the differential rotation zone was taken to be ~ 1/2 
the distance from the base of the convective zone to the stellar core. In addition, AQ and 
Q were assumed to be constant and independent of the magnetic field. Thus, the dynamo 
lasted indefinitely, and sustained a toroidal field of ~ 5 x 10'' G at the interface. The 
arbitrarily long lifetime was essential because in order for the large-scale field to drive a 
self-coUimated outfiow, the dynamo must last until the end of the AGB phase after the 
star has radiatively bled most of its envelope material. To alleviate the assumption of 
steady AQ and ft and to study the backreaction of field growth on the shear, we apply our 
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dynamical dynamo model using parameters in Blackman et al. (2001). 

We fix the thickness of the differential rotation zone (L = 4.6 x 10^*^ cm), the rotational 
speed at the interface (Qq = 5 x 10~^ rad s""^), and the shear profile across the Q-layer. 
We then lower the value of /3<^. This lengthens the dynamo lifetime and determines the 
relative fraction of energy deposited into magnetic or heat sinks. Fig. 2 shows the isolated 
interface dynamo for two different values of f3^. In both cases the depth of the differential 
rotation zone is fixed and the shear energy in that zone is available for field amplification. 
The toroidal and poloidal fields are out of phase and thus generate an oscillatory Poynting 
fiux. Maximum toroidal field strengths of ~ 5 x lO'^ G are obtained, while the poloidal field 
is significantly lower, with values of i?p ~ 2 x 10^ G at the interface. In both cases, the 
time integrated Poynting flux and turbulent dissipation generate ~ 10~^ times the binding 
energy of the envelope. Therefore, the current parameters can not produce a dynamically 
influential magnetic outflow or signiflcant heating from turbulent dissipation in the shear 
zone. If an isolated interface dynamo were to be viable, a deeper shear zone would be 
required. 

A higher c<^ allows the poloidal component of the held to diffuse deeper into the 
toroidal zone, thus extracting more differential rotation energy. Hence, the depth of the 
shear layer is determined by the value of c^. If we fix L as the distance from the base of the 
convective zone to the core, then the corresponding fields diffuse all the way to the stellar 
core for Pvp > 10~®. This unbinds the envelope at the beginning of the AGB phase. In 
addition, we can constrain how far the fields would have to diffuse in order to unbind the 
AGB envelope. If the poloidal field diffuses to ~ 5 x 10^ cm within the stellar interior, then 
Maq = 1-7 X 10^1 g and Afi = 1.5 x 10"=^ rad s-^ This corresponds to E^n ^ 1-5 x 10^^ 
ergs which is comparable to the binding energy of the envelope. However, in these cases, the 
dynamo would blow off the envelope prematurely and such a circumstance would contradict 
observations of AGB hfetimes ~ 10^ yrs. If an isolated shell dynamo is to be consistent 
with observations, the dynamo must be sustained during the end of the AGB phase. 



4.2. Isolated Dynamo With Reseeding AQ 

Convection might reinforce differential rotation analogous to what occurs in the sun. 
Thermal energy and a negative entropy gradient drive convective turbulence. In the 
Kolmogorov framework, energy from the large-scale to the dissipative scale cascades at a 
rate given by 

(-) 
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where I is a length scale, vi the corresponding velocity and D is independent of scale (Shu 
(1992)). For our PNe progenitor, we use the large-scale convective velocity, v = vi = 10^, 
cm/s corresponding to the approximate thickness of the convective zone, I = Li. We 
envision adding a resupply term to Eq. (10), in which a fraction, /, of the turbulent energy 
cascade rate resupplies shear. To determine this term, we equate (15) with the time rate of 
change of kinetic energy in the convective zone and arrive at 



where Afl^ is the shear that would be resupplicd by the turbulent cascade. By adding the 
right side of Eq. (16) to Eq. (10), we can investigate the full range of convective resupply 
scenarios from / = (no convective resupply) to / = 1 (maximum convective resupply). 
In addition, we fix the following values: AQq = 1.5 x 10~^ rad/s, Qq = 5 x 10~® rad/s, 
L = 4.6 X 10^° cm, Cp = IQ-^ and Q = 5.0. 

We consider two sub-scenarios for the convective resupply dynamo: Q dynamically 
evolving and Q constant. In the first case, Poynting fiux is allowed to spin down the 
envelope at the interface. This requires magnetic buoyancy, which appears in Eq. (4) as 
terms proportional to uB/L. The left graph in Fig. 3 presents the result for dynamical 
Q,. Rotation is drained by Poynting flux and we take the most extreme best case of 100 % 
(/ = 1) of the turbulent energy cascade rate resupplying differential rotation energy. This 
does establish a constant Afi. However, the bulk of the Poynting flux is drained in a short 
burst (< 10 yrs). It is unfeasible to generate the requisite energy required to produce a 
magnetically dominated explosion in this case. 

For the second case (Q constant), we consider the possibility that the field is stored in 
the interface layer until the corresponding aggregate Poynting flux is able to blow up the 
envelope through a magnetic "spring" effect. If the field is trapped, Poynting flux does not 
emerge from the layer and thus does not spin down the envelope. Even though the dynamo 
equations for the two cases are mathematically identical, in the case of steady rotation 
terms proportional to uB/L can be interpreted as diffusion rather than buoyant loss. The 
right plot in Fig. 3 shows a constant Q with / = (no convective resupply). Constant 
rotation results in a decay of the differential rotation energy. However, since Q is constant, 
the a-effect is non-zero and thus is continually pumping poloidal fleld into the shear layer. 
The achieved poloidal fleld strength is neghgible and energetically insigniflcant. It cannot 
blow off the envelope. 

Finally we consider the case in which AQ is resupplied and Q, is constant. Fig. 4 
demonstrates that convective resupply coupled with a constant Q results in a sustained 
Poynting flux. The Poynting flux of ~ 5 x 10^^ erg/s sustained over a period of 10^ yrs is 




(16) 
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enough to overcome the binding energy of the envelope and produce a magnetically driven 
envelope expulsion. For = 10~^, 0.1% of the energy cascade rate has been used to 
resupply the differential rotation and marks the approximate minimum threshold fraction 
required to blow off the envelope. If > 5 x 10^^, then the rate of heat produced from 
turbulent dissipation is greater than the Poynting flux, resulting in thermally induced 
envelope expulsion. 

Although it is not known whether a mechanism for resupplying differential rotation 
operates in AGB stars (e.g. Riidiger & Hollcrbach (2004)), we have demonstrated that such 
an effect may facilitate a dynamo driven envelope expulsion, provided both Q and AQ are 
sustained. 



4.3. Common Envelope Dynamo 

In a close binary system, Roche lobe overflow can result in both companions immersed 
inside a common envelope (Paczynski (1976), Ibcn & Livio (1993)). A drag force due to 
the velocity difference between companion and envelope, induces in-spiral of the secondary. 
As a result, orbital energy and angular momentum are transfered from companion to 
common envelope. For low-mass secondaries, the in-spiral time can be extremely fast (< 1 
yr) supplying the requisite angular momentum on time scales much shorter then the AGB 
hfetime (Nordhaus & Blackman (2006)). As the companion traverses the envelope, the 
transfer of orbital energy alone may be enough to unbind it. However, if the companion 
cannot supply the necessary orbital energy, in-spiral continues until the secondary is tidally 
shredded into a disk. In addition, angular momentum transfer spins up the envelope 
resulting in a differentially rotating stellar interior. Assuming that the AGB star is initially 
stationary, we can calculate the corresponding rotation profile generated by in-spiral of a 
companion. The change in orbital energy, in virial equilibrium, of the secondary is given as 

AB„,(r) = ^Ml!^ - (17) 

where Mt is the total mass of the star, m2 the secondary mass, M the enclosed mass at 
position r, and is the stellar radius. A fraction, aAEorb of the orbital energy released by 
the companion is available for mass ejection. If aAEorb > Ebind, where Ebind is the envelope 
binding energy, then the secondary has expelled the envelope. In order to calculate the 
rotation profile of the AGB star, we use the gravitational potential energy released by the 
secondary to spin-up spherical shells as follows: 



aAEorb = h^l (n) , 



(18) 
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where Ig is the moment of inertia of a thin shell, Qg the angular velocity of the shell and 
the outer radius of the shell. Solving for the rotation profile of a shell yields 



where Mg is the shell mass. 

We study the limit in which the orbital eccentricity is negligible, resulting in the 
companion exhibiting Kcplerian motion at all radii (Pollard (1979)). Fig. 5 shows envelope 
rotation profiles for companions with masses ranging from 0.01 to 0.05 Mq. Also shown 
is the Keplerian velocity, vk, and the sound speed, Cg inside the envelope. Higher mass 
companions supply enough orbital energy to spin-up the envelope above its Keplerian value 
at a given radius. The rotational energy will then redistribute via outward mass transfer 
until Keplerian rotation is re-established. This effect can be seen in the envelope rotation 
profile generated by a 0.05 Mq secondary. It also occurs for larger mass companions, 
however we do not investigate those as we would expect similar results. The dashed-dotted 
vertical line shows the approximate radius at which the companion is tidally shredded and 
spin-up of the envelope ceases. Below this depth the companion would form a disk and 
the mechanism of angular momentum transfer would be different. In reality the rotation 
profile should be solved for self-consistently, but this becomes particularly important as 
soon as the rotational energy exceeds that associated with the sound speed. We have not 
incorporated this non-linear effect explicitly, and hence our approach of redistributing the 
excess rotational energy in the inner regions in Fig. 5 is crude. 

As can be seen from Fig. 5, a range of companion masses can produce varying amounts 
of rotation and differential rotation. To investigate how this shear energy is deposited, 
we apply the rotation profiles in Fig. 5 along with the parameters of the stellar model to 
our interface dynamo. For a 0.02 Mq brown dwarf secondary, toroidal field strengths of 
~ 1 X 10^ G arc obtained (see Fig. 6). For Pvp = 10^"^, the decay of the shear energy 
and toroidal field are long (~ 25 yrs) and occur over several thousand cycle periods. 
Therefore, in Fig. 6, the solid line represents the envelope of the dynamo while the smaller 
insets show a "zoomed- in" region to demonstrate the oscillatory nature of the fields. 
For these parameters, the heat generated from turbulent dissipation is greater then the 
time-integrated Poynting fiux, thus we identify this as a thermally induced model. 

In Fig. 7, we present a model in which the time-integrated Poynting fiux is large 
enough to unbind the stellar envelope. In this case, the companion is a 0.05 Mq brown 
dwarf and spins-up the envelope until it is shredded into a disk. For this model, Prp — 10~^ 
with the corresponding Poynting fiux decaying in ~ 100 yrs. Peak toroidal and poloidal 
field strengths are comparable to results from Fig. 6, however the lower Prp results in 
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less differential rotation energy being converted into heat. Therefore, we identify this as a 
magnetically dominated model. 

Both magnetically dominated and thermally driven models can be produced for a range 
of companion masses and diffusion coefficients. The resultant outflows for the two cases are 
expected to be qualitatively different. For interface-dynamo-driven winds, the launching 

and shaping of the outflow could occur close to the core. Such an outflow is expected to be 
collimated, predominately poloidal and aligned with the central rotation axis. However, if 
heat is the primary driver mediating the transition from progenitor to PPNe, the resulting 
outflow is probably quasi-spherical. Bipolar, magnetically collimated PNe could be the 
result of our low Pr^, common envelope magnetically dominated models. 

5. Conclusions 

Extraction of rotational energy is likely fundamental to the formation of multipolar 
PPNe and PNe. Magnetic dynamos can play an intermediary role in facilitating the 
extraction of rotational energy. Here, we have studied a large-scale, dynamical interface 
dynamo operating in the envelope of a 3.0 Mq AGB star. The back reaction of fleld 
amplification on the shear is included as are the drain of rotation and differential rotation via 
both turbulent dissipation (thermally induced) and Poynting flux (magnetically induced). 
Two different dynamos are studied: (i) that of an isolated star and (ii) that of a common 
envelope system in which the secondary is a low-mass companion (< 0.05 Mq). 

For the isolated case, we flnd that only when two conditions are met can the single 
star dynamo drive PPNe. First, AQ must be re-seeded. This may occur by analogy to the 
A-effect in the sun (Riidiger & HoUerbach (2004)). Secondly, Q must be constant. This 
implies that the field is stored in the shear layer until the end of the AGB phase. When 
these two stringent conditions are met, a small fraction of the energy cascade rate can 
provide the necessary shear energy to sustain the interface dynamo (in some cases as little 
as 0.1%). Not only is the dynamo maintained, but the time-integrated Poynting flux is 
large enough to overcome the envelope binding energy. 

Whether or not isolated star dynamos can produce a PPNe, a binary interaction can do 
so more robustly for a wide range of cases. Common envelope evolution is advantageous in 
several ways. Energy and angular momentum are supplied very quickly (< 1 yr) and often 
in less than or equal to a dynamo cycle period, allowing the dynamo to operate once the 
secondary has completed its in-spiral. For our common envelope dynamo model, a range 
of companion masses can easily supply enough differential rotation energy to power either 
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a dynamo driven jet or a thermally driven outflow. A magnetically dominated explosion 
likely produces a coUimated, poloidal outflow while that of a thermally induced explosion is 
expected to be more spherical. 

We have highhghted some of the basic key issues of isolated and common envelope 
dynamos, however, more detailed research is needed in both areas. The viability of 
anisotropic convection reseeding differential rotation must be determined. If convection 
cannot reseed shear, then we are faced with the proposition that binary interactions are 
required to produce axisymmctric PNc. In the CE case, the complex interaction between 
companion and envelope, multi-dimensional aspects of the dynamo, realistic rotation 
profiles of isolated stars and the inclusion of a wider array of secondary masses are just 
a few of the many problems which warrant future work. Constraints on the turbulent 
diffusion coefficient as a function of radius also need to be determined. The physics involved 
in transitioning from a dynamo to a fully active jet must be understood. 
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Fig. 1. — A meridonial slice of the dynamo geometry. The left figure shows the global 
geometry of the AGB star. The right figure is a close-up view of the dashed region on the 
left. The a-effect is driven by convection and occurs in layer of thickness Li above the 
differential rotation zone. The poloidal component of the field is pumped downwards into 
the differential zone, where it is wrapped torodially due to the fi-effect. 
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Fig. 2. — The differential rotation energy is allowed to drain through field amplification 
and turbulent dissipation. In this figure, A;^ = 5 x 10~^^ cm~^, Cp — 0.01, Q — 5.0. We 

define, e[PF,dis] = lo E[PF^dis] (f) dt' and label M (epp) and T (e^is) on the top right plot to 
distinguish between the thermal and magnetic contributions to the binding energy. For the 
left figure, = 10^^ while the right has = 10^^. Peak field strengths are a factor of 
~ 5 — 10 less then those obtained in (Blackman et al. (2001)). Differential rotation energy 
is drained in < 20 yrs. Lowering results in the differential rotation energy draining at a 
slower rate, allowing the field to sustain for longer periods of time (~ 40 — 50 yrs). However, 
peak field strengths remain the same. 
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Fig. 3. — Results for reseeding differential rotation through convection. In the left figure, 
/ = 1 corresponding to maximum convective resupply. Rotation is drained through Poynting 
flux but cannot sustain a dynamically important dynamo. In the figure on the right, / = 
(no resupply of Afl). The rotation rate is fixed, corresponding to a buildup of Poynting flux 
in the interface layer. 
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Fig. 4. — Convective resupply results in a steady-state differential rotation profile. For the 
left column, the envelope of the poloidal, toroidal and Poynting flux is plotted. The Poynting 
flux is sustained at ~ 5 x 10^^ erg/s. The sustained Poynting flux supplies enough energy to 
unbind the envelope of our 3 Mq model at the end of the AGB phase (~ 10^ yrs). In this 
flgure, C0 = 10~^ and / = 10~'^ implying that only ~ 0.1% of the cascade energy must be 
converted into differential rotation energy to supply the requisite Poynting flux. This model 
predicts a magnetically dominated explosion. 



Fig. 5. — Rotation profiles generated from the transfer of angular momentum from 
companion to envelope in our AGB star. In both figures, the solid curves represent the 
resulting profiles for companions of masses 0.05 (top), 0.02 and 0.01 (bottom) Mq. For the 
left figure, the dashed curve represents the Keplerian velocity while the dashed-dotted curve 
is the sound speed. The 0.05 Mq companion initially spins up the envelope such that the 
inner region is rotating faster then the Keplerian velocity. Mass redistribution ensues and 
transfers matter outward until the rotation profile drops below Keplerian. The right figure 
presents the angular velocity corresponding to the left figure. The dash-dot vertical line is 
the approximate radius at which the companion is tidally shredded. The large-dash vertical 
line is the boundary of the shear layer in Blackman et al. (2001) while the small dash fine 
is the base of the convection zone. These profiles assume that a — 0.3. 
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Fig. 6. — Interface dynamo resulting from the in-spiral of a 0.02 Mq brown dwarf in the 
interior of our model AGB star. The differential rotation zone extends from the base of 
the convection zone to the radius at which the secondary is tidally shredded (Nordhaus 
& Blackman (2006)). In this model, Pm = 10~^ and Q ^ 5, Qo = 2.3 x 10"^ rad/s, 
AQo = 2.5 X 10~^ rad/s and 5/L — 1. In the left column, the envelope of the Poynting flux 
(top), toroidal field (middle) and poloidal field (bottom) are drawn with a solid line. The 
insets represent the time evolution from to 0.2 yrs. The vertical scale of the insets are the 
same as the corresponding larger figure. 
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Fig. 7. — Interface dynamo resulting from the in-spiral of a 0.05 Mq brown dwarf. In this 
model Pm = 10-^ Q = 5, 6/ L = 1, = 5 x 10"^ rad/s and Af^o = 2.5 x lO^^ rad/s. The 
insets represent the time evolution from to 0.2 years. The vertical scale of the insets are 
the same as the corresponding larger figure. 



